Observation of the critical regime near Anderson localization of light 



in 
o 
o 

(N 
> 

o 

(N 



S3 



c3 



i 

C3 
O 

o 



(N 
> 

oo 



in 
o 

c3 



C3 
O 

o 



X 



Martin Storzer, Peter Gross, Christof M. Aegerter, and Georg Maret 
Fachbereich Physik, University of Konstanz, Germany 
(February 2, 2008) 

Diffusive transport is among the most common phenomena in nature [1]. However, as predicted 
by Anderson [2], diffusion may break down due to interference. This transition from diffusive 
transport to localization of waves should occur for any type of classical or quantum wave in any 
media as long as the wavelength becomes comparable to the transport mean free path £* [3]. The 
signatures of localization and those of absorption, or bound states, can however be similar, such 
that an unequivocal proof of the existence of wave localization in disordered bulk materials is still 
lacking. Here we present measurements of time resolved non-classical diffusion of visible light in 
strongly scattering samples, which cannot be explained by absorption, sample geometry or reduction 
in transport velocity. Deviations from classical diffusion increase strongly with decreasing £* as 
expected for a phase transition. This constitutes an experimental realization of the critical regime 
in the approach to Anderson localization. 



Diffusive transport controls many phenomena in bi- 
ology, chemistry, physics and engineering, where undi- 
rected transport is characterized by a linear increase of 
the mean square displacement (r 2 ) with time [1]. As 
a consequence, transmission of particles through a slab 
of length L is proportional to £* / L known in the con- 
text of electrons in metals as Ohm's law. However, when 
the wave-nature of the diffusing particles is taken into 
account, constructive interference of waves propagating 
on reciprocal multiple scattering paths may lead to a 
breakdown of diffusion and the particles are "trapped" 
on closed multiple scattering loops. This means that the 
probability of returning to the starting point of such a 
closed loop is increased twofold due to the fact that the 
path has the same length in both counter-propagating di- 
rections. Thus diffusive transport away from this starting 
point is reduced. Anderson first predicted this destruc- 
tion of diffusion in 1958 [2] to explain the metal-insulator 
transition. If the scattering power of the medium is high, 
i.e. the transport mean free path I* is small, the density 
of closed loops increases which leads to a reduced trans- 
port inside the material. This can be described by a 
rescaling of the diffusion coefficient [4] . For strong local- 
ization the rescaled diffusion coefficient becomes zero at 
finite length scales, such that the spread of (r 2 ) comes 
to an end at a length scale corresponding to the localiza- 
tion length. The transition from a diffusive to a localized 
state should occur when the wavelength A becomes com- 
parable to I*, as quantified by the Ioffe-Regel criterion 
kl* ~ 1 where k = 2ir/\ is the wavenumber [3]. 

While localization of electrons in metals was widely 
studied [5], an unequivocal demonstration of the pure 
interference effect as predicted by Anderson turned out 
to be difficult in this electrons bound in deep 

minima of a random potential are virtually impossible to 
distinguish from those in closed loops. However, photons 
in disordered materials revealed an excellent model sys- 
tem to study localization effects [6,7]: In contrast to elec- 



trons, they do not interact with each other, nor can they 
be bound in a random potential. Hence all deviations 
from classical diffusion have to be due to interference, 
when absorption is taken into account properly. More- 
over, in order to create a strongly scattering medium, ma- 
terials having a large optical refractive index combined 
with very small absorption are available for visible light. 

Thus there have been many investigations in the past 
studying diffusive transport of visible [8,9] and infrared 
light [10]. For visible light, where reported values of kl* 
go down to 3.2 [8], measurements have mainly focused 
on static transmission experiments where a decay faster 
than oc X/L was interpreted as evidence for the onset 
of localization [12]. It was noted however that absorp- 
tion, which is always present in such samples, can lead 
to a similar decay [13,14]. Furthermore, recent investiga- 
tions of time resolved transmission of these samples with 
similar values of k£* show no discernible deviation from 
purely diffusive behavior [9] . Localization of photons has 
also been studied by use of microwave scattering in dis- 
ordered media [14-16], where indications of a decreasing 
diffusion coefficient have been found [16]. However it 
has to be noted that these investigations were carried 
out in tubes of extent ~ 2 x 2 x 20 in units of £* [17], 
which constitute a relatively small quasi one-dimensional 
(ID) system. As localization is always present in ID and 
2D [4], theoretical descriptions of the quasi ID geometry 
[17,18] indicate that observations of localization effects 
are possible far away from k£* ~ 1. Similar theories in 
3D [19] find a negligible effect in the same regime. Thus, 
an unequivocal proof of the transition to localization in 
three dimensions is still lacking. 

Here we present measurements of time resolved pho- 
ton transport through bulk powders of Ti02 with typical 
sample dimensions of ~ 10 5 x 10 5 x 10 4 in units of I* . 
These measurements have the advantage that they allow 
a direct determination of the photon path length depen- 
dence of the diffusion coefficient, where absorption and 
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localization lead to different functional forms of the tail 
(exponential vs. non-exponential) of the path length dis- 
tribution [18]. Furthermore we independently measure 
the turbidity of the samples using coherent backscatter- 
ing [20,21]. This allows a systematic study of the ap- 
proach to the localization transition as a function of tur- 
bidity (hi*)- 1 . 

Our samples consist of ground TiC>2 particles in its ru- 
tile structure with a refractive index of 2.8. The resulting 
grains of the different samples have various average par- 
ticle sizes with average diameters ranging from 220 nm 
to 550 nm with a polydispersity of ~ 25%. Those par- 
ticles are commercially available as pigments for white 
paint from e.g. DuPont or Aldrich. In order to mini- 
mize I* these powders are compressed to obtain filling 
fractions of <f> ~ 0.4. The scattering properties are deter- 
mined from coherent backscattering [20,21]. The angular 
width of this enhancement to the incoherent background 
is inversely proportional to kl* [22,23]. For the highly 
scattering samples studied here, the width becomes very 
broad, which is why the distribution has to be measured 
to very wide angles for a good determination of the inco- 
herent background. Our setup was custom designed for 
this problem (see methods) [24]. 

Fig.l shows the backscattering angular distribution for 
two different samples, SI and S3, corresponding to aver- 
age particle diameters of 250 nm and 550 nm respec- 
tively. Using the value of the averaged refractive index 
from numerical simulations based on the energy density 
coherent potential approximation [25], we determine the 
reflectivity of the surface and are hence able to correct 
for the overestimation of kl* due to internal reflections 
[26]. Thus the values of kl* are determined from the 
full widths at half maximum of the curves with a correc- 
tion for the reflectivity of the surface. This gives values 
ranging from kl* = 6.3 for S3 to kl* = 2.5 for SI at a 
wavelength of 590 nm, close to the Ioffc-Regel criterion. 

In order to measure the time resolved transmission we 
use a single photon counting method, where the time de- 
lay of a picosecond light pulse transmitted through the 
sample is measured. From a histogram of time of flights 
of many such pulses, the path length distribution inside 
the sample is obtained directly [9,27,28]. Our setup con- 
sists of a dye laser working at a wavelength of 590 nm 
with a pulse width of ~20 ps. The time of flight his- 
tograms have been deconvoluted with the pulse shape of 
the laser system in order to recover the pure path length 
distributions (see methods). Fig. 2 shows the path length 
distributions in the samples characterized in Fig. 1 com- 
pared to another sample (S2) with ki* = 4.3. As can be 
seen in Fig. 2A, classical diffusion theory including ab- 
sorption [29,30] fits the data from sample S3 very well. 
In contrast, the path length distribution of SI (Fig. 2C) 
shows marked deviations from diffusive behaviour. In 
particular, we observe a non-exponential decay at long 
times, with photons staying inside the sample longer than 



expected from a purely diffusive process. In Fig. 2B, the 
data from sample S2 show small deviations from the dif- 
fusion picture. 

Diffusion in an extended slab of thickness L, which in 
our experiments ranges from 1.3 mm to 2.5 mm, leads to 
an exponential decay of the time dependent transmission 

( 7T 2 -D(t) | e V 

at long times as I(t) oc e V L ™ a J [29,30]. Here 
l a is the absorption length, which in our samples ranges 
from 0.3 m to 2.6 m, c is the speed of light in vacuum, n 
is the effective refractive index of the medium, and D the 
diffusion coefficient. Note that the values for l a are more 
than a factor of 10 6 larger than typical values of I* in 
our samples. Any non-exponential decay of I{t) at long 
times thus indicates a temporally varying diffusion coef- 
ficient D{t). Using the properties of I(t) at long times, 
we can obtain a direct measurement of the diffusion coef- 
ficient up to a constant given by l a . This is achieved by 
taking the negative time derivative of the logarithm of 
the intensity divided by (-K 2 D{t = 0)/L 2 + c/nl a ), where 
n 2 D(t = 0)/L 2 is the inverse of the time t maxi at which 
most photons leave the sample [16]. This is shown in 
Fig. 3 for samples SI, S2, and S3. For long times the 
classical sample S3 approaches the constant value of 1, 
whereas the curve for SI shows a strongly decreasing dif- 
fusion coefficient and S2 reveals intermediate behaviour. 

To quantify the non-classicality of the diffusive light 
transport through the different samples we determined 
the average of the ratio of the measured data to the dif- 
fusion fits. This was done systematically for all powders 
over the time interval from t max to 3t max . This devia- 
tion is given as a function of kt* in Fig. 4. The figure 
clearly shows that the deviations increase strongly when 
approaching the Ioffe-Regel criterion, as expected for a 
phase transition. 

The data we have presented show clear deviations from 
diffusive transport through these highly scattering sam- 
ples. The non-exponential decay of the path length distri- 
bution indicates a renormalized value of the diffusion co- 
efficient at long times as was predicted by scaling theory 
[4,7]. Moreover, the deviations from diffusive behaviour 
scale with the value of kl* consistent with the Ioffc- 
Regel criterion [3] and the approach to a phase transition. 
These deviations cannot be explained by absorption as 
this only leads to an additional exponential decrease but 
cannot introduce a non-exponential path length distri- 
bution. There is also no systematic dependence of the 
values of the absorption length of the various samples on 
kl* in contrast to the monotonic dependence of the devi- 
ations from classical diffusion. A change in the transport 
velocity due to resonant scattering [31] can also be ruled 
out, since values of the diffusion coefficient from the clas- 
sical fits, such as that in Fig. 2C, are consistent with 
a value of the transport velocity of c/n. Furthermore 
such a change is not expected to be dependent on the 
length of the specific path and hence would not lead to a 
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non-exponential decay as the data show. Finally, sample 
inhomogeneities, such as a stratification, cannot account 
for the effect, as the path length distributions are found 
to be independent of the direction of illumination. We 
thus conclude that these observations constitute direct 
evidence for the slowing down of photon diffusion due to 
the approach to the Anderson localization transition and 
hence for the existence of the transition to strong local- 
ization of photons in three dimensions. 



METHODS 
Coherent backscattering 

The diffuse reflected intensity is not constant but de- 
pends on the backscattering angle 9 as cos(#), such that 
the difference of the backscattering signal from this func- 
tional shape has to be determined. Our setup consists 
of 256 photo-sensitive diodes attached to an arc with a 
diameter of 1.2 m in order to get sufficient angular reso- 
lution over a range of -70° < 9 < +70° [24]. Here, the 
resolution is 1° for \9\ > 20°, 0.7° for 10° < \6\ < 20° and 
0.14° for \9\ < 10°. In addition, the central part of the 
backscattering cone, —3° < 9 < +3° was measured sepa- 
rately using a beam-splitter and a charged coupled device 
camera to a resolution of 0.02°. The measurements are 
done using circularly polarized light in order to reduce 
the influence of singly scattered light. 

Time of flight measurements 

Our setup consists of a dye laser working at a wave- 
length of 590 nm with a pulse width of ~20 ps. The 
Rhodamin6G dye laser (Coherent C699) is pumped by an 
Ar + laser (Coherent Inova400) and was modified by APE 
with a mode locker and a cavity dumper to deliver pi- 
cosecond pulses. In order to recover the pure path length 
distributions, the time of flight histograms have to be de- 
convoluted with the pulse shape of the laser system. This 
is because after pulses, albeit strongly suppressed, and in- 
discriminate noise may lead to disturbances in the time 
of flight measurements. For this purpose, we measured 
pulse shape, including the noise level, in the absence of 
a sample. This zero-pulse was then deconvoluted with 
the time of flight data in Fourier space to directly give 
the path length distribution for a supposed delta-peaked 
pulse as it is calculated theoretically. 

The observed non-exponential decay is not due to a 
distribution of I* values in the sample due to e.g. strati- 
fication. We have confirmed this by measuring the time 
of flight distributions with sample SI flipped, such that 
the direction of stratification would be reversed. These 



measurements show the same non-exponential decay as 
Fig. 2C. 
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FIG. 1. Measurements of coherent backscattering for two 
different samples. Open symbols: sample SI with an average 
particle diameter of 250 nm which yields kt = 2.5. Closed 
symbols: sample S3 with an average diameter of 550 nm and 
kt — 6.3. All measurements were done with circularly polar- 
ized light at a wavelength 2jr/k = 590 nm. 
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FIG. 2. Path length distributions from samples SI, S2, and 
S3. The experimental results are compared to diffusion theory 
including absorption. In A, one can see that the data from 
S3 (L = 2.5 mm, D = 22 m 2 /s, £ a = 2600 mm) closely follow 
the diffusion fit, showing an exponential decay at long times. 
Part C in contrast shows strong deviations from the diffusion 
fit for SI (L = 1.48 mm, D — 15 m 2 /s, l a = 340 mm), with a 
clearly non-exponential decay at long times. These deviations 
can be explained by a time dependent diffusion coefficient in 
the sample. An intermediate case is shown in part B from 
sample S2 (L — 1.51 mm, D = 13 m 2 /s, t a = 380 mm), with 
a value of kl* = 4.3, where small deviations from the classical 
behaviour can be observed. 
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FIG. 3. Taking the negative derivative of the logarithm of 
the measured data one obtains a measure of the time depen- 
dence of the diffusion coefficient out of the time of flight data. 
Sample S3 has a constant diffusion coefficient at long times 
in agreement with diffusion theory indicated by the solid line. 
SI and S2 in contrast show a decay illustrating the decrease of 
the diffusion coefficient of light with time. This demonstrates 
the existence of localized modes. 
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FIG. 4. Measurements of the deviations from diffusion on 
several samples with different values of kl*. The deviation 
is quantified by taking the ratio of the experimental path 
length distribution, I(t), to the classical fit, If it (t), and av- 
eraging it over the time interval from t max to 3t max . The 
results were plotted over the k£* value determined from co- 
herent backscattering (see text). Samples that have a kl* 
nearer to the Ioffe-Regel criterion show a bigger deviation 
from diffusive behaviour. 



5 



